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We study the existence of chimera states in pulse-coupled networks of bursting Hindmarsh-Rose 
neurons with nonlocal, global and local (nearest neighbor) couplings. Through a linear stability 
analysis, we discuss the behavior of stability function in the incoherent (i.e. disorder), coherent, 
chimera and multi-chimera states. Surprisingly, we find that chimera and multi-chimera states 
occur even using local nearest neighbor interaction in a network of identical bursting neurons alone. 
This is in contrast with the existence of chimera states in populations of nonlocally or globally 
coupled oscillators. A chemical synaptic coupling function is used which plays a key role in the 
emergence of chimera states in bursting neurons. Existence of chimera, multi-chimera, coherent and 
disordered states are confirmed by means of the recently introduced statistical measures and mean 
phase velocity. 


I. INTRODUCTION 

Synchronization in networks of neuronal systems has 
been an active research area due to its important role in 
coding and information processing in biological systems 
and brain. Bursting neurons are characterized by alter¬ 
nates of activity of neurons between a quiescent state 
and fast repetitive spiking on a slow time scale. There 
are many processes to produce bursting behavior in cou¬ 
pled oscillators [1]. Further, different types of synchrony 
occur in coupled bursting neurons which include spike 
synchronization, bursting synchronization, complete syn¬ 
chronization and anti-phase bursting synchronization [2]. 
In general, burst synchronization occurs at lower values 
of coupling strength and for complete synchronization, 
which involves both spike and burst synchronizations, 
one requires a higher value of coupling strength. Two 
different forms of couplings are mainly used in coupled 
bursting neurons depending on whether the synapse is 
chemical Q or electrical Q. In the first case, the synap¬ 
tic coupling is often approximated by a static sigmoidal 
nonlinear input-output function with a threshold and sat¬ 
uration. In the second case, the electrical coupling is a 
linear function and is directly dependent on the difference 
between the membrane potentials. Complete synchro¬ 
nization always occurs in globally coupled bursting neu¬ 
rons under suitable conditions on the coupling strength 
when the synaptic coupling is electrical. Using such a 
global synaptic coupling, identical oscillators are found 
to be either synchronized or oscillate incoherently, but 
they never exist simultaneously. Interestingly, we observe 
the coexistence of coherent and disordered (incoherent) 
states under nonlocal, global as well as local (nearest 
neighbor) synaptic couplings in a network of bursting 
neurons. Such a coexisting state was originally named as 
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chimera state by Abrams and Strogatz [5] in the context 
of nonlocally coupled phase oscillators. 

Chimera states in identical coupled oscillators have 
also been an active area of extensive research in recent 
years in the field of biology, physics and social sciences 
[3, 16]. Initially the chimera states have been shown to 
emerge when a network of identical oscillators are cou¬ 
pled nonlocally, that is the coupling strength decays with 
distance between the coupled oscillators. Chimera states 
are interesting because they occur even when the oscil¬ 
lators are identical and coupled symmetrically. There 
are different types of chimera states such as amplitude- 
mediated [7], breathing |j|, clustered [9j, spiral-wave 0, 
etc. types. Chimera states have also been experimentally 
observed in chemical [Ill, electronic 112], electrochem¬ 
ical HI and mechanical 0 systems. Initially it was 
assumed that chimera states only exist in phase oscilla¬ 
tors using non-local coupling configuration. But recently 
it has been observed that chimera states also occur in 
systems exhibiting limit cycles and in chaotic dynamical 
systems 0. Very recently, chimera states have also been 
observed in globally coupled oscillators [l6|. Chimera 
states observed in real-world systems 0 , where various 
dynamical behaviors are involved. For example, in the 
case of Parkinson’s disease due to lose or damage cells 
in the brain, synchronized activity is absent in certain 
regions of the brain 0. In the case of epileptic seizures, 
specific regions of the brain are highly .synchronized and 
the others part are not synchronized [l9| . 

One of the most challenging and inspiring problems in 
this area is to identify the existence of chimera states 
in neurobiological systems. During the last decade the 
mechanism behind mutual synchronization and phase 
synchronization of chaotic bursts in neural ensembles has 
been explained in detail sin eh. Recently, Kalitzin et 
al.p2j have observed collective dynamics of coupled neu¬ 
ronal oscillators which have multiple oscillatory states. 
In this paper we report a surprising find of chimera and 
multi-chimera states in networks of chaotically burst¬ 
ing Hindmarsh-Rose oscillators [H] under different cou- 
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plings, including local ones. Chimera states emerge in an 
ensemble of bursting neurons with at least three types 
of couplings, namely nonlocal, global and local (near¬ 
est neighbor) interactions. Recently, Hizanidis et al. 
0 investigated the existence of chimera states in three- 
dimensional Hindmarsh-Rose neuron model using a non¬ 
local electric type of coupling. Previously, chimera states 
were observed by Sakaguchi [25] in coupled Hodgkin- 
Huxley neural oscillators with excitatory and inhibitory 
couplings where nonlocal coupling is essential for the ap¬ 
pearance of these states. Very recently, Omelchenko et 
al. [26] have discussed the robustness of chimera states 
in nonlocally coupled FitzHugh-Nagumo oscillators with 
respect to perturbations of the frequencies of the in¬ 
dividual oscillators and the structural transformations 
of the network topology. Belykh et al. [2j| discussed 
the influence of coupling strength and network topology 
on synchronization in pulse-coupled network of bursting 
Hindmarsh-Rose neurons, where no chimera like behavior 
has been reported. 

As noted above, chimera states have been observed in 
Hindmarsh-Rose neurons with electric type of nonlocal 
coupling [24|, where the range of coupling and coupling 
strength play key roles. But the existence of chimera 
state in Hindmarsh-Rose neurons with chemical synap¬ 
tic coupling under local, global and nonlocal interactions 
has not been reported earlier. In our studies a chemical 
synaptic coupling is used which plays a key role in the 
emergence of chimera states in bursting neurons for lo¬ 
cally and globally coupled neurons. We also analytically 
derive the condition for mutual synchronization in glob¬ 
ally pulse-coupled network of bursting neurons. In the 
case of global and local (nearest neighbor) coupling, the 
behavior of the stability function in incoherent, chimera 
and multi-chimera states are also discussed. In our anal¬ 
ysis we have used suitable statistical measures recently 
introduced by Gopal et al. [27] as well as the notion of 
mean phase velocity proposed by Omelchenko et al. [28] 
to confirm the chimera and multi-chimera states. 

The subsequent part of this paper is organized as fol¬ 
lows. Sec. II is devoted to a brief presentation of 
Hindmarsh-Rose model for a bursting neuron. In Sec. 
Ill, we numerically investigate the existence of chimera 
states using nonlocal coupling. In sec. IV, an analyti¬ 
cal study of synchronization in globally coupled neurons 
is reported. The behavior of the stability function in 
chimera states is also discussed. Existence of chimera 
state in local (nearest neighbor) coupling is described in 
Sec. V. Sec. VI provides a discussion of our results. 

II. HINDMARSH-ROSE NEURON MODEL 

The Hindmarsh-Rose neuron model [23], which is a 
well known system for its chaotic behavior and differ¬ 
ent types of bursting, in its original form is expressed as 
follows: 

x = y + ax 2 — x 3 — z V /, 


y = 1 - dx 2 - y, (1) 

i = c(b(x — xo) — z), 

where the variable x represents the membrane potential, 
and the variables y and 2 are the transport of ions across 
the membrane through the fast and slow channels, re¬ 
spectively. The fast variable y represents the rate of 
change of the fast (e.g. sodium) current. The slow vari¬ 
able z controls the speed of variation of the slow (e.g. 
potassium) current. This speed is controlled by the small 
parameter c. The parameter I describes an external cur¬ 
rent that enters the neuron and xo is a control parameter 
delaying and advancing the activation of the slow cur¬ 
rent in the modeled neuron. For the sake of simplicity, 
after the linear transformation / parametric redefinition 
j29| x —y x, y —y 1 — y •> z —y IT I z ^ d —y a V cq ^ ^ 

— 1 — I — bx 0 , Eq.(l) can be written in the form 

x = ax 2 — x 3 — y — z, 

V = (a + a)x 2 - y, (2) 

i = c(bx — z T e). 

This transformed model (2) is a phenomenological model 
that can exhibit all common dynamical features found 
in a number of biophysical modelling studies of burst¬ 
ing. We consider c a small positive parameter so that 
z(t) varies much slower than x(t) and y(t). Square-wave 
bursting has been observed for the following set of pa¬ 
rameter values: a = 2.8, a = 1.6, c = 0.001,5 = 9 and 
e = 5 [20]. The system (2) is monostable, that is coexis¬ 
tence of stable equilibrium point and limit cycle has not 
been observed for this set of parameter values. 

III. NON-LOCAL INTERACTION 

Now we consider a network of identical Hindmarsh- 
Rose neurons with non-local interaction as follows: 

^ j=i+P 

Xi = ax1 -xj - Vi- Zi + —(v s - Xi) ^2 CijT(xj), 

j=i—p 

Vi = (a + a)x • - y u 

(3) 

Zi = c(bxi~Zi+e), i = 1, 2, •••, N 

where N is the total number of elements in the network, 
p is the number of coupled nearest neighbors in each di¬ 
rection on a ring so that the coupling radius r = and k 
is the synaptic coupling strength. The connectivity ma¬ 
trix C = ( Cij) nxn is such that q? = 1 if the i th neuron 
is connected to the j th neuron, otherwise Cij = 0 and 
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ca = 0. The synapse is excitatory for the reversal poten¬ 
tial v s = 2 > Xi(t ) for all times t and Xi(t). The synaptic 
coupling function T(x) is modeled by the sigmoidal non¬ 
linear input-output function as 


1 _|_ g—A(x—© s ) 5 

where A determines the slope of the function and © s is 
the synaptic threshold. This often-used coupling form 
has been called the fast threshold modulation [3j|. We 
choose the threshold © s = —0.25 so that every spike in 
the isolated neuron burst can reach the threshold. We 
fixed the value of A = 10 throughout the work. From a 
physicist’s perspective, equation (3) represents a network 
of TV identical pulse-couped oscillators with nonlocal at¬ 
tractive interaction. But from a neuroscientist’s point 
of view, such a network corresponds to an interaction 
between 7V-nonlocally coupled excitatory neurons with 
direct excitatory synapses [30|. 



FIG. 1. (Color online) Left panels show snapshots of ampli¬ 
tudes for (a) disordered state at k = 0.3, (b) chimera state 
at k = 0.85 and (c) coherent state at k = 1.4. Right panels 
(e), (f) and (g) show the corresponding mean phase velocity 
cUi,i = 1,2, • • -,7V, for each neuron. The coupling radius is 
fixed at r — 0.3, where TV — 200. 

Now, we consider numerically the existence of chimera 
states and transitions of chimera and multi-chimera 
states as well as disordered and coherent states in the 
network of nonlocally coupled Hindmarsh-Rose oscilla¬ 
tors. The network (3) and their variants in the follow¬ 
ing sections are integrated using the fifth order Runge- 
Kutta-Fehlberg integration algorithm scheme with inte¬ 
gration step length At = 0.01. The initial conditions 
for (3) are chosen as follows: Xio = 0.01(7 — y),^o = 
0.02(7 - f),* i0 = 0.03(7 - f) for 7 = 1,2,- ■ -,f 
and for rest of the oscillators = 0.1 — i),Vio = 
0.12(y — i),Zio = 0.21(y — 7) with added small ran¬ 
dom fluctuations. In Fig. 1, we fix the coupling radius 
as r = 0.3 and vary the synaptic coupling strength k. 


In the left panel, the snapshots of the amplitude (mem¬ 
brane potential) for disordered, chimera, and coherent 
states for k = 0.3,0.85, and 1.4, respectively, are shown. 
To confirm chimera or coherent states in nonlocally cou¬ 
pled neurons, first we calculate the mean phase velocity 
recently proposed by Omelchenko et al. [28| of each neu¬ 
ron as uji = 2irMi/AT, i = 1,2,3, • • -,7V, where Mi is the 
number of bursts of the 7th neuron during a sufficiently 
long time interval AT. Figure 1 (right panel) shows the 
mean phase velocity of each neuron corresponding to in¬ 
coherent, chimera and coherent states for the synaptic 
coupling strength k = 0.3,0.85 and 1.4 respectively. To 
calculate the mean phase velocity, the time interval is 
taken over 4 x 10 5 time units after an initial transient of 
1 x 10 5 units. 

To clearly distinguish the disordered, chimera, multi¬ 
chimera and coherent states, we also use the recently 
introduced statistical measures by Gopal et al. [22j us¬ 
ing the time series of the network. For this purpose we 
will calculate the strength of incoherence (SI) and dis¬ 
continuity measure (DM) from a local standard devia¬ 
tion analysis. To calculate these statistical measures, we 
first introduce a transformations = Xi — a^+i, W 2 ,i = 
m - Vi+hW^i = Zi - Zi+i, 7 = 1,2, •• -,7V [see Ref. 
[27j for more details]. Different synchronization states 
in the network can be quantified by using the standard 
deviations given by 





N 

Ei»w - wp 


( 5 ) 


where l = 1,2,3; z = 1,2, • • -,7V, and (wi) = 

~k ^ZiLi an d (• • -) t denotes the average over time. 

For coherent states the values of o\ s are zero while they 
take non-zero values for both disordered and chimera 
states. To distinguish chimera and disordered states, we 
divide the number of oscillators into M (even) bins of 
equal length n = N/M. Then we introduce local stan¬ 
dard deviation which is defined as 


< n(m ) 



1 

n 


mn 

E Ki - WP 

j=n(m— 1) + 1 


(6) 


m = 1, 2, • • *, M. The above quantity 07 ( 771 ) can be calcu¬ 
lated for every successive n number of oscillators. Then 
the strength of incoherence (SI) is defined as 


SI = 1 — 


A/m= 1 


Sm 


M 


s m = Q(S - cri(m )), 


( 7 ) 


where ©(•) is the Heaviside step function and S is a pre¬ 
defined threshold. Consequently, the values of SI=1 or 
SI=0 or 0 < SI < 1 represent disordered, coherent and 
chimera or multi-chimera states respectively. Again in 
order to distinguish chimera and multi-chimera states we 
also introduce the so called discontinuity measure (DM) 
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[27| which is defined as 

DM = 1 Y 1 — with s M+ 1 = Sl . (8) 

For chimera state the value of DM = 1 and for multi¬ 
chimera state the value of DM is a positive integer greater 


than “1”. 
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FIG. 2. (Color online) Strength of incoherence SI and dis¬ 
continuity measure DM are plotted against synaptic coupling 
strength k for two values of coupling radius r. (a, b) and 
(c, d) correspond to r — 0.3 and r — 0.4, respectively. Here 
N = 200, M = 40 and S = 0.05. 


Fig. 2 shows the variation of the strength of incoher¬ 
ence and discontinuity measure for fixed value of cou¬ 
pling radius r and different values of synaptic coupling 
strength k. To start with we choose N = 200 and the 
total number of bins to be M = 40 which we find to be 
optimal. For r = 0.3, the variation of SI (Fig. 2(a)) 
and DM (Fig. 2(b)) are shown for different values of k. 
To calculate standard deviation 07 in Eq. (5) and local 
standard deviation 07 ( 771 ) in Eq. (6), the time average is 
taken over t = 4 x 10 5 time units after an initial tran¬ 
sient of 10 5 units. The existence of disordered, chimera, 
multi-chimera and coherent states are represented by SI 
and DM. As mentioned above, (SI, DM) = (1, 0) repre¬ 
sents a disordered state, while (SI, DM)=(0, 0) repre¬ 
sents a coherent state. Further 0 < SI < 1, DM=1 and 
0 < SI < 1, 2 < DM < M/2 represent chimera and multi¬ 
chimera states, respectively [23]. For r = 0.3, at lower 
values of fc, we observe that all the neurons in the net¬ 
works are in a disordered state, represented by the region 
I ={k : 0 < k < 0.72}. With an increase in the value of k 
beyond k = 0.72 we observe chimera and multi-chimera 
states in the region ll={k : 0.72 < k < 1.24}. Chimera 
states represented by the value of DM = 1 and multi¬ 
chimera states with DM = 2, 3,4, 5, 6 are shown in Fig. 
2(b). In this state, all the neurons, namely coherent and 
disordered groups are in bursting type. With a further 


increase in the coupling strength we observe that all the 
neurons of the network are to be in coherent states rep¬ 
resented by the region III={/c : 1.24 < k < 2.0}. For 
k > 2.0, all the neurons are always in a coherent state. 
Similarly, the variations of SI and DM for different values 
of the synaptic coupling strength k are shown in Fig. 2(c) 
and Fig. 2(d), respectively, for r = 0.4. From Fig. 2, we 
observe that the extent of the region II for chimera and 
multi-chimera states gets decreased with increasing value 
of the coupling radius r. Fig. 3 shows the two-parameter 
(r, k) phase diagram of coherent, chimera / multi-chimera 
and incoherent states. The range of synaptic coupling 
strength k for chimera or multi-chimera states is large 
for small values of coupling radius r, that is the network 
is almost locally coupled. But this region decreases for 
higher values of r near to 0.5, that is when the system 
is globally coupled. Finally we have confirmed the above 
results for larger sizes of the network, namely N = 300 
and N = 500 neurons to make sure that our results hold 
good for larger networks as well. 



FIG. 3. (Color online) Two parameter (r, k) phase diagram for 
N = 200 nonlocally coupled network of identical Hindmarsh- 
Rose oscillators. Strength of incoherence is used as a mea¬ 
sure for incoherence, coherence and chimera / multi-chimera 
states. Red region is for coherent, green is for chimera or 
multi-chimera and gray is for incoherent states. 


IV. GLOBAL INTERACTION 


We can convert the non-local interaction into a global 
interaction by taking the number of nearest neighbors 
as p = where the number of nodes N in the net¬ 

work is an odd number. Then the network of identical 
Hindmarsh-Rose oscillators with global coupling is as fol- 
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lows: 

k * 

Xi = ax--Xi~yi-Zi+ - (v s -Xi)22c ij T(x j ), 

3 =i 

Vi = (a + a)x 2 - y h 

(9) 

Zi = c(b Xi~ Zi + e ), i = 1,2, • • •, TV 

where C = ( Cij) nxn is the connectivity matrix such that 
Qj = 1 if i 7 ^ j and cu = 0. 

In order to understand the existence of mutual 
synchronization of the above globally coupled net¬ 
work, we carry out an analytical investigation based 
on a linear stability analysis, closely following the 
stability analysis for a synchronized manifold by 
Belykh et al. 0 m. For the synchronized state 
Xi(t) = x(t),yi(t) = y(t),Zi(t) = z(t) for all i, and so the 
system (9) becomes 

X = ax 2 -x 3 -y-z+ j^- - (v s -x)T(x), 

y = (a + a)x 2 -y , 

(10) 

i = c(b x—z+e), 

where R is the row sum of the connectivity matrix C. 
Introducing the differences between the neural oscillator 
coordinates, 

€ij =x j -x i ,rn j —yj-yutij =Zj-z h iJ = 1,2, ••-,7V, 

(ii) 

the linearized stability equations for the transverse 
perturbation of the synchronization manifold are given 

by 


iij — (2 ax 3x ryj Qj jy _ ^ 

k , 

■b jy _ ^ ( v s — ^)F x ( x ){Rfiij + / X C jh^jh — Cihfiih)), 

h= 1 

rjij = 2 (a+a)x Zij-rjij, 

(12) 

Cij = c(b 

where x(t) is the synchronous solution defined by the 
system (10) and T a ,(x) is the partial derivative of T(x) 


with respect to x. The term ^2h=i( c Mjh ~ Cih£ih) is 
the same as in the case of linear coupling |2(J, [32| . The 
derivatives are calculated at£ = 0,77 = 0,£ = 0. Then 
the stability equations become 

£ = (2ax-3x 2 )£-r]-(-n(x)£, 

f] = 2 (a+a)x 77, 

(13) 

C = c(b£- C), 

where fi(x) = ~ 7T=i( v » ~ x ) T x( x )( R + ^ 2 ), A 2 

is the largest real part of the eigenvalues of the coupling 
matrix M = C — RI. It is well known that the matrix 
M has one zero eigenvalue Ai and all the other eigenval¬ 
ues have non-positive real parts [20]. If the coupling is 
mutual then the coupling matrix M is symmetric and all 
the eigenvalues are real. For simplicity, let us suppose 
that the largest eigenvalue A 2 of the coupling matrix M 
is simple. Equation (13) is then an analog of the Master 
Stability equation [1321 ]. 


(a) (b) 
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FIG. 4. (Color online) Snapshots of a system of globally 
coupled Hindmarsh-Rose neurons for different values of the 
synaptic coupling strength k in terms of the variables Xi 
(black color) and the transformed variables wi,i — Xi — Xi+i 
(red color): (a) incoherent state, k — 1.0, (b) chimera state 
(with two synchronized/desynchronized groups), which can 
be reindexed into single group each, k — 1.2, (c) chimera 
state (with single synchronized group), k = 1.28, (d) coherent 
state, k = 1.3. The inset figures are the corresponding time 
series (blue color). The number of oscillators is N = 301. 

Numerical results of the global synaptically coupled 
Hindmarsh-Rose neurons are presented in Fig. 4 for 
N = 301. The initial conditions are chosen as follows: 
X i0 = 0.01(i - ^i), y i0 = 0.02(i - ^fi), z i0 = 0.03 (i - 
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for i = 1,2,-**, and x* 0 = 0.1(^ J: - 2/*o = 

0.12(^^ -i),^ 0 = 0.21 (A^I -7) for 7 = ^=^- + 1, •••,7V 
with added small random fluctuations. The snapshot of 
the state variables x^ and the new transformed variables 
wi,i = Xi — Xi +1 are shown by black and red color dotted 
points, respectively. At a lower value of the synaptic cou¬ 
pling strength k = 1.0, all the neurons are in a disordered 
state (Fig. 4(a)) and their behavior are all of square-wave 
bursting in nature. At a higher value of k = 1.2, a typi¬ 
cal pattern of the chimera state with two groups of syn¬ 
chronized/desynchronized oscillators, which of course can 
always be brought into separate single group of synchro¬ 
nized/desynchronized oscillators by appropriate reindex¬ 
ing of the oscillators (due to the global nature of the 
coupling), is observed between the neurons (Fig. 4(b)). 
The behavior of a neuron in the coherent group and a 
neuron in the disordered group have the same time se¬ 
ries form, that is coexistence of square-wave and plateau 
bursting is observed in all the neurons. The typical time 
series of x^ (blue color line) is shown in the inset of Fig. 
4(b). The values of the transformed variables for the 
coherent groups are near zero whereas for the disordered 
neurons they are randomly distributed in [-2, 2]. With 
a further increase in the value of the synaptic coupling 
strength to k = 1.28, chimera state with single coher¬ 
ent and incoherent group each is observed in which the 
left group of neurons are in coherence where the trans¬ 
formed variables take values close to zero and the right 
one is an incoherent one (Fig. 4(c)). All the coherent 
and disordered neurons exhibit a similar behavior, that 
is a mixture of square-wave and plateau bursting states 
shown in the inset of Fig. 4(c) by blue color line. At a 
further higher value of the synaptic coupling strength at 
k = 1.3, all the neurons are found to be in a coherent or 
completely synchronized state (Fig. 4(d)). In this case 
all the neurons are in plateau bursting states. In globally 
coupled neurons, the region of chimera states are very 
small compared to non-local coupling, as already shown 
in Fig. 3 for r = 0.5 and N = 200. 

To confirm chimera states in globally coupled neu¬ 
rons, we calculate the mean phase velocity uji(i = 
1,2,3,- • •, TV) of each neuron. Figure 5 shows the mean 
phase velocity of each neuron corresponding to incoher¬ 
ent, chimera (with two synchronized / desynchronized 
groups), chimera (with single synchronized group) and 
coherent states. We have numerically integrated the dy¬ 
namical equations (9) of the globally coupled neurons and 
the initial conditions used are the same as above. The 
time interval is taken over t = 5 x 10 5 time units after an 
initial transient of 10 5 units. We obtained similar results 
for even larger number of neurons, namely N=401 and 
501, which confirm our above conclusions. 

The stability function for mutual synchronization is 
n(x) = - x )T' x (x)(R + V), with R = 

TV—1 and A 2 = —TV. It is strongly dependent on the mem¬ 
brane potential x(t). The synaptic coupling strength k is 
the upper bound of O(x) for x(t) > 0 S and rapidly de¬ 
creases to zero if x(t) < 0 S . The global coupling yields a 
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FIG. 5. Mean phase velocities ca of globally coupled neurons 
corresponding to Fig. 4: (a) disordered state at k — 1.0, 
(b) chimera state (with two synchronized / desynchronized 
groups) at k = 1.2, (c) chimera state (with single synchronized 
group) at k = 1.28, and (d) coherent state at k = 1.3. Other 
parameters are same as in Fig. 4. 
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FIG. 6. (Color online) Variation of stability function Q(x) 
(green color), nullcline z — f(x) (red color) and phase tra¬ 
jectory projected in x — z plane (blue color) are shown for 
(a) disordered state at k — 1.0, (b) multi-chimera state at 
k — 1.2, (c) chimera state at k = 1.28, and (d) coherent state 
at k — 1.3. 

nullcline of x as 2 : = f(x) = —ax 2 — x 3 — x)T(x). 

The behavior of the stability function 0(x) and null- 
cline 2 = f(x) at incoherent, chimera (with two syn¬ 
chronized / desynchronized groups), chimera (with sin¬ 
gle synchronized group) and coherent states are shown 
in Figs. 6(a - d). It is important to note that at lower 
values of synaptic coupling strength, square-wave burst¬ 
ing is observed which turns into plateau bursting states 
for higher values of k when coherent states occur in the 
globally coupled excitatory network (9). This happens 
when the synaptic coupling strength is large enough to 
change the square-wave to plateau bursting through a 
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homoclinic bifurcation. By changing the synaptic cou¬ 
pling strength the system (9) undergoes transition from 
disordered to coherent states corresponding to the tran¬ 
sition from square-wave to plateau bursting state via the 
disappearance of the homoclinic bifurcation. During the 
transition regime chimera states are observed. In the 
chimera states a combined bursting of square-wave and 
plateau is observed. 


V. LOCAL INTERACTION 

Next we consider a network of identical Hindmarsh- 
Rose oscillators with local (nearest neighbor ) interac¬ 
tion, that is p = 1 in the non-local interaction (3), as 
follows: 

k 

Xi = axf-Xi-yi-z i +-(v s -x i )(r(x i - 1 )+r(x i+1 )), 


m = {a+a)xi~yi, 


(14) 


Zi = c(b Xi~Zi+e), i = 1,2, •••,7V, 

with periodic boundary conditions (xq = xn,xn+i = 
xi). As before, here k is the synaptic coupling strength. 

A sequence of interesting behaviors occur numerically 
by changing the synaptic coupling strength k. We used 
the same initial conditions as was used in the case of non- 
locally coupled systems earlier in Sec. III. Fig. 7 shows 
the snapshots of the collective behavior of neurons for dif¬ 
ferent values of the coupling strength k. We observe that 
at lower values of the coupling strength, all the neurons 
are in disordered states. This is shown in Fig. 7(a) for 
k = 0.4. With increasing value of k and beyond a critical 
value, we observe multi-chimera state in Fig. 7(b) for 
k = 1.2. With further increase of &, this multi-chimera 
state is transformed into a chimera state in Fig. 7(c) for 
k = 1.36. For a further increase in the coupling strength 
to k = 3.6 all the neurons are in a coherent state [Fig. 
7(d)]. 

To further clarify the collective behavior, we plot the 
spatio-temporal nature of all the nodes xi(t) in the right 
panel of Fig. 7. Color coding of Fig. 7(e), 7(f), 7(g) 
and 7(h) clearly shows the presence of disordered, multi¬ 
chimera, chimera and coherent states for four distinct 
values of the synaptic coupling strength k = 0.4,1.2,1.36, 
and 3.6, respectively. Snapshots of the locally cou¬ 
pled neurons (14) for different values of the synaptic 
coupling strength in terms of the new state variables 
wij = Xi — Xi +1 are shown in Fig. 8. From this figure it 
is seen that initially at a lower value of the synaptic cou¬ 
pling strength, k = 0.4, all the neurons are in a disordered 
state and the values of are randomly distributed over 
[-2, 2] in Fig. 8(a). On increasing the synaptic coupling 
strength to k m 1.2 [Fig. 8(b)] and k = 1.36 [Fig. 8(c)], 
the values of for coherent group of neurons are near 



FIG. 7. (Color online) Snapshots of amplitude (left panel) 
for (a) disordered state at k — 0.4, (b) multi-chimera state at 
k — 1.2, (c) chimera state at k = 1.36, and (d) coherent state 
at k — 3.6. Figs, (e)-(h) show spatio-temporal color coded 
maps for (a)-(d), respectively. 



i i 


FIG. 8. (Color online) Snapshots of the transformed variables 
wi,i corresponding to Fig. 7. (a) disordered state, k m 0.4, 
(b) multi-chimera state, k — 1.2, (c) chimera state, k = 1.36 
and (d) coherent state, k = 3.6. 


zero but the disordered group of neurons are randomly 
distributed over [-2, 2]. On further increasing the cou¬ 
pling strength to k = 3.6, [Fig. 8(d)] all the neurons are 
in a coherent state and the values of w\^ are near zero 
for all times. 

To confirm disordered, multi-chimera, chimera and co¬ 
herent states, we have first calculated the mean phase 
velocity uji for each neuron which are shown in Fig. 9. 
The mean phase velocity uji for each neuron in disor¬ 
dered, multi-chimera, chimera and coherent states for 
k = 0.4,1.2,1.36 and 3.6 are shown in Figs. 9(a)-9(d) 
respectively. The values of mean phase velocity uji in 
Fig. 9(d) for each neuron is near to zero for the coher¬ 
ent state (k = 3.6). In the coherent state the behavior of 
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all the neurons is close to the steady state of the locally 
coupled system x = ax 2 — x 3 — y — z + k(v s — x)T(x),y 
(a + a)x 2 — y, z = c(b x — z + e). 



FIG. 9. Values of mean phase velocity uj% for each neuron 
corresponding to Fig. 7: (a) disordered state, k = 0.4, (b) 
multi-chimera state, k — 1.2, (c) chimera state, k m 1.36 and 
(d) coherent state, k = 3.6. Other value of parameters are 
same as Fig. 7. 


In order to distinguish multi-chimera and chimera 
states clearly we use the strength of incoherence and dis¬ 
continuity measure for different values of synaptic cou¬ 
pling strength k. In Figs. 10(a) and 10(b) we demonstrate 
the behaviors of the strength of incoherence and discon¬ 
tinuity measure as a function of the synaptic coupling 
strength k which clearly confirms the presence of chimera 
and multi-chimera states in the case of local interaction 
also. This appears to be quite surprising in the sense that 
chimera states have normally been identified only in the 
case of nonlocally or globally coupled arrays. Existence of 
chimera states using local coupling has been observed for 
larger number of neurons also, namely N = (300,500). 




FIG. 10. (Color online) Variation of (a) strength of inco¬ 
herence (SI) and (b) discontinuity measure (DM) for differ¬ 
ent values of synaptic coupling strength k for locally cou¬ 
pled Hindmarsh-Rose neurons. Here N = 200, M — 40 and 
S s= 0.05. The time average t — 5000 is considered after an 
initial transient of 10 5 time units to calculate cri and cri(m) in 
Eqs. (5) and (6). 


We also note that for local interaction in Eq. (3) 
Cij = 1 for either j = i + 1 or j = i — 1 so that 
(14) is valid in this case. For all other values of j, 
= 0 and cu = 0. The value of eigenvalue A2 for a 
ring of 2iF-nearest neighbor mutually coupled neurons is 
A 2 = —4 YltLi sin 2 ^ In the case of local (nearest 

neighbor) interaction (Eq. (14)), K = 1 and so the value 
of A2 is A2 = —4 sin 2 (-^) with N = 200. 

Now considering the master stability function |32| for 
the local synaptically coupled network (14), the stabil¬ 
ity function is Q(x) = RkT(x) — k(v s — x)T x (x)(2 + A 2 ) 
with R = 2. For local excitatory coupling the nullcline 
of x is z = f(x) = —a x 2 — x 3 + ^f-(v s — x)T(x). 



FIG. 11. (Color online) The function £}(#) (green color line), 
nullcline f{pc) (red color line) and corresponding projected 
phase trajectories in x — z plane (blue color line), (a) Inco¬ 
herent bursting for k = 0.4 and (b) multi-chimera state for 
k = 1.2. The behaviour of the neuron (black color line) at 
k — 0.4 and k = 1.2 are shown in inset figures of (a) and (b) 
respectively. 


The function Q(x) corresponds to the contribution of 
the synaptic coupling for the stability of the synchroniza¬ 
tion manifold. In fact, in the case of global coupling the 
stability function Q(x) varies between zero and synaptic 
coupling strength k for all x(t). On the other hand for 
local coupling Q(x) <0 for all values of the membrane 
potential x(t) < ©s and, while it has an upper bound k 
for x(t) > O s 0 , |3l|. The first term of Q(x) contains 
the Heaviside-type synaptic function T(x) and it becomes 
significant when x(t) > 0 S . The second term of Q,(x) is 
decisive for the values of x near the threshold © s as the 
derivative T x has a rapidly decaying tail. Practically the 
cells are uncoupled when the value of x(t) is below the 
threshold 0 S . 

At a lower value of the synaptic coupling strength 
k = 0.4, Q(x) has a negative drop near x = 0 S = —0.25 
(Fig. 11(a)) which signifies the incoherent state of the 
bursting neuron. The corresponding time series of the 
bursting neuron is shown in the inset of Fig. 11(a). With 
a further increase of coupling at k = 1.2, the value of 
Q(x) is negative and has a sharp drop near x = 0 S . The 
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neurons are in a multi-chimera state. The coherent and 
disordered neurons are in a plateau bursting state in Fig. 
11(b) which is in contrast to the chimera states stud¬ 
ied earlier in the literature in chaotic systems, where the 
coherent states are in periodic states or remain close to 
steady states and disordered states are in chaotic states 
j34j . But in a chimera or multi-chimera state in the 
present case all the coherent as well as the disordered 
neurons are in plateau bursting states. 

VI. CONCLUSION 

In conclusion, we have analyzed the role of the chem¬ 
ical synaptic coupling function in inducing chimera and 
multi-chimera states of bursting neurons using non-local, 
global and local (nearest neighbor) interaction of neu¬ 
rons. Surprisingly, we find that chimera or multi-chimera 
states occur even in the presence of local interaction 
alone whereas previous studies |J [6] of chimera or multi¬ 
chimera states exist using either non-local or global in¬ 
teraction only. Interestingly, we identified three types of 
chimera states using different types of coupling configura¬ 
tions. In the first type, in the case of non-local coupling, 
in chimera or multi-chimera states, coherent and disor¬ 
dered neurons are all of bursting type. In the second 
type, we find that in the chimera or multi-chimera state 
using global synaptic coupling both the groups, namely 
the coherent or disordered groups, are in a combination 
of square-wave and plateau bursting states. That is for 
some times they are square-wave bursting in nature and 
for other times they are in a plateau bursting state (as 


shown in Fig. 4(b) or 4(c)). In the third type of chimera 
state using local synaptic coupling, we identified that 
both the groups in chimera / multi-chimera states are in 
plateau bursting states (Fig. 11(b)). Using suitable sta¬ 
tistical measures disordered, multi-chimera, chimera or 
coherent states are confirmed. The existence of chimera 
and coherent states are also confirmed by using mean 
phase velocities [28]. To conclude, we wish to point 
out that this work promises to identity the existence of 
chimera states in various types of coupling topologies in 
bursting cells. Some types of chimera states might be 
helpful for information processing in the case of neuro¬ 
logical diseases. Existence of chimera states in two di¬ 
mensional grid of oscillators is more realistic in neurobi¬ 
ology. In our manuscript we observed chimera state using 
nonlocal, global and local type of interactions. Since in¬ 
teraction in two dimensional grid of oscillators are also 
of these types, we believe that chimera states can also 
emerge in these cases too. We are presently exploring 
this phenomenon and hope to report the results in the 
near future. 

Finally, it will be also of interest to extend our work to 
identify chimera or multi-chimera states in other systems 
using nonlinear local coupling. Moreover, the question of 
introducing time-delay in the nonlinear local coupling for 
the existence of chimera states [35j is also an important 
task for future studies. 
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